Dynamical model for the full stretching curve of DNA 
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We present a phenomenological dynamical model able to describe the stretching features of a 
length vs applied force DNA curve. As concerning the chain, the model grounds on the discrete 
worm-like chain model with the elastic modifications, which properly describes the elongation fea- 
tures at low and intermediate forces. The dynamics is developed under a double well potential with a 
linear term, which, at high forces, accounts for the narrow transition present in the DNA elongation 
(overstretching). A quite good agreement between simulation and experiment is obtained. 

PACS numbers: 87.15.-v, 36.20.-r, 87.18.Tt, 83.10.Rs, 05.40.-a 



Introduction. The stretching curve of DNA has been 
the subject of many studies, both theoretical and ex- 
perimental. The first experimental study of the DNA 
stretching was performed by Bustamante and co-workers 
[l[ , where the DNA length was measured as a function of 
the force applied to the chain. From the theoretical point 
of view, the small and intermediate force range of the 
DNA stretching curve can be understood and depicted 
by means of the wormlike chain (WLC) model Q, where 
the DNA is modeled as a flexible beam. For simulation 
purposes, the WLC model can be discretized as a chain 
of beads connected by sticks with the presence of an elas- 
tic bending (discrete WLC) and improves the more naive 
freely joined chain model Q). This WLC model (and the 
discrete version) fits with excellent precision the experi- 
mental curve, making use of the correction introduced by 
Odijk [|| that takes into account the extensibility of the 
chain when the sticks are replaced by harmonic springs. 

The most intriguing feature of the stretching of DNA 
is presented in a subsequent experiment of Bustamante 
and co-workers [H, @ , where a sudden elongation of the 
chain is registered when a large force, around 70pN, is 
applied. 

The presence of this tension-induced overstretching 
transition reveals the existence of two structurally differ- 
ent DNA states. The first state, at low applied forces, is 
the so called B-DNA conformation, where the base pairs 
arc hclicoidally packed in their native state. The other 
state, present at high forces, was first called the S state, 
where the base pair distance is higher by about 70% than 
that in the state B. The nature of the overstretched DNA 
state remains controversial, and it seems that many ef- 
fects could be responsible for the sudden DNA elongation 
at large forces. 

Molecular dynamics simulations revealed that the 
DNA transition could correspond to a change of the he- 
licoidal DNA shape to another one with a shorter helix 
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radius and a large tilt of the base pairs [6J. This result 
seems also be confirmed by early spectroscopic studies 

0- 

On the other hand, the suggestion that the enlarge- 
ment is accompanied by a rotation of the bases giving 
rise to a ladderlike structure in the S state of DNA, has 
also been made J8|. Motivated by the experimental data 
of Strick et al. [9[ , a possible linear coupling of the twist- 
stretch variable was proposed [HI [ll|, an d in this last 
paper an estimation of the DNA twist stiffness was also 
performed. More accurate measurements suggest that a 
not completely laddcrlike form of the DNA occurs in the 
elongation from the B- to the overstretched DNA state. 
In fact a remanent helicity of the DNA persists in its over- 
stretched state of about one turn every 37.5 base pairs, 
almost four times smaller than the natural helicity of one 
turn per 10.5 base pairs [l2|. Moreover, by avoiding the 
writhing of the DNA double helix, a broadening in the 
overstretching transition is observed. 

Various statistical-mechanics models have been imple- 
mented in order to mathematically describe the narrow 
overstretchin g tr ansition, which take into account the 
WLC model [La - [l4j . with an increasing number of pa- 
rameters in order to take into account the twist con- 
straints dud. 

The existence of the S-DNA state as a hybridized state 
of DNA has been the object of deep and still controver- 
sial investigation, mainly arising because of the presence 
of an asymmetric hysteresis in the pulling and relax- 
ing DNA elongation curve, which appeared to indicate 
a force- induced DNA melting fl6l - [Tc| . Further investiga- 
tions, both numerical and experimental, were published 
up until very recently, some of which confirmed fl9l - [22| , 
and some did not confirm j23M25| this hypothesis. 

In this work, we present a dynamical Langevin model 
where a polymer is pulled under an asymmetric free en- 
ergy potential, with the purpose of describing the narrow 
DNA overstretched transition at high forces in a Landau- 
Ginzburg landscape. This model proposes a simple phe- 
nomenological way to describe the DNA dynamics, and 
its results are in good agreement with the most relevant 
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FIG. 1: Scheme of linear chain in 3D space. In (a) the chain 
is free to move in the thermal bath, while in (b) is pulled from 
each end with a force F v . 
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FIG. 2: Potential for different values of the force F p acting 
on the chain. The forces increase from the upper to the lower 
curves. 



experimental outcomes. The model we present doesn't 
enter in the microscopic details of the DNA elongated 
state, if melting is or is not responsible for it. The de- 
scription presented is valid regardless of the specific mi- 
croscopic state of the DNA. In this sense the model is 
general and can be included as an ingredient into other 
more complex processes like DNA translocation driven 
by molecular and nanotechnological motors. 

The chain model. We model the DNA molecule by a 
three-dimensional (3D) polymeric chain of N dimension- 
less monomers connected by nonharmonic springs (JT|) . In 
spite of the usual harmonic interaction |26f . our clastic 
potential energy is given by 
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where k e is the elastic parameter, is the position of the 
i-th particle, = |d;| = |rj + i — is the distance be- 
tween the monomers i + 1 and i, and Iq and l\ are minima 
of the potential representing approximately the equilib- 
rium distance between adjacent monomers, at weak and 
strong forces respectively. The potential is plotted in 
Fig.[5J The parameter fc; is chosen in such a way that the 
middle of the force transition corresponds to equal prob- 
abilities to cross the elongation potential barrier from 
cither left to right (enlargement) or right to left (con- 
striction). 

The model takes into account the bending energy with 
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FIG. 3: DNA length vs Force. Parameters: t s = 5 x 10 6 , 
N = 100, k B T = 4.1, k b = 10, k e = 1190, and h = 89 (units 
given in the text). 
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where kb is the bending elastic constant, and 6i is the 
angle between the link dj+i and the link d^. 

The dynamics of the chain is given by the overdamped 
equation of motion 



Ti = -ViV e} (di) - V i y ben (e i ) - y/2k B T£(t), (3) 

where £ (t) represents the thermal contribution as a Gaus- 
sian uncorrelated noise, and the time t is scaled with the 
damping 7 as t —> t/"f. Time units are then given in 
S7 -1 . The nabla operator is defined as Vj = d/dxi'i + 



d/dyij + d/dzik. Two forces F p and 



-F p act on the 



chain on the first and last monomer respectively in order 
to stretch the polymer. The figures show the length of 
the DNA polymer as a function of the applied force F p . 

Results. We performed various computer simulations 
with a Runge-Kutta stochastic algorithm for a long fixed 
simulation time t s = 5 x 10 6 time units (t.u.), and 
dt = 0.01. We set the units of energy and length to 
be E u = 4.1 pN nm (thermal energy at room tem- 
perature) and l u = 5.3 nm, respectively With these 
units the parameters used in simulations are ksT = E u , 
k e = U90EJI*., l = lu, h = l.7l u , h = 10EJ U , 
and 9 — 0. Thus the persistence length of the chain 
is £ = kt/ksT = I0l u = 53nm, which gives the correct 
WLC model fit from the experimental data. The param- 
eter k e is determined from the extended WLC model in 
order to fit the DNA elongation data in region II, and 
analogously for the bending parameter kb- 

After a thermalization time of t s = 1 x 10 6 t.u., the 
average of the lengths at all subsequent times has been 
evaluated with respect to the contour length d c of the 
polymer, with d c — (N — 1)Z , and the values are plot- 
ted in Fig. [3J where a very good agreement between 
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the computer simulation (brown □) and the experiments 
(black x) is shown 

The inset of the figure shows good agreement also at 
low force values. The full line represents the prediction 
of the WLC model in region I and II, which follows very 
well both the simulations points and the experimental 
ones. Figure H] shows the density distribution of the link 
length close to the transition. It is there evident that 
a gradual increase of the number of enlarged links (i.e., 
di »s 1.7) is produced by increasing the applied pull force 
F p . 

The transition between the regimes II (B-DNA) and IV 
(overstertched) arises in the very short force range. We 
model this feature using the Ginzburg-Landau (GL) ef- 
fective potential given in Eq. ([IJ . This approach has also 
been used to describe the B-DNA to Z-DNA transition 
[2^| . In the GL phase transition theory an asymmetrical 
contribution to the free energy gives rise to a discontin- 
uous (first-order) phase transition with the temperature 
T as control parameter [28| . Analogously, in the present 
case, we obtain a sharp transition with the DNA elonga- 
tion d x as order parameter, and the force F p as control, 
with the difference with respect to the standard behav- 
ior that the intensity of thermal fluctuations is here an 
external parameter, whose presence causes a relatively 
smooth transition between the two states. 

A dynamical fingerprint of discontinuous transitions 
is the appearance of hysteretic behavior due to a slow 
convergence to equilibrium states. This behavior has 
been clearly observed in out-of-equilibrium experiments 
on the overstretching transition [23|, [U, (29|. The hys- 
teresis curves are consequence of force ramps applied on 
a time scale shorter than the DNA relaxation time. Using 
waiting times of the order of 1 min the equilibrium state 
is recovered and no hysteresis is observed [L2j . Figure [5] 
shows how our model mimics this behavior. In fact, the 
left panel of the figure presents a typical hysteresis curve 
obtained by decreasing the simulation time, i.e., not al- 
lowing the system to reach equilibrium. For a simulation 
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FIG. 4: Probability density of the length of the links di for 
different applied forces along the overstretching transition. 
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FIG. 5: Left panel: transition region with different applica- 
tion times of the force, t B = 3 x 10 J , and t s = 10 6 , shorter 
than the relaxation time of the chain. Right panel: the same 
for different values of the potential parameter ki , ki increasing 
from left to right curves. 

time of t s = 10 6 t.u., we already have a different path 
for increasing forces than for decreasing ones (indicated 
in the figures by the black arrows). For a shorter simula- 
tion time t s = 3 x 10 5 t.u., the hysteresis area is bigger, 
while for t s — 5 x 10 6 t.u. (squares in the figure), no 
hysteresis is observed. This behavior expresses a typical 
dynamical property of the model not visible in a static 
equilibrium analysis. 

Finally, we discuss the effect of the ki parameter on 
the transition. The right panel of Fig. [5] shows how this 
parameter is responsible in the model for the force value 
where the transition occurs. This value can be related 
to those parameters whose magnitude affects the over- 
stretching transition: the value of pH [IB] or the salt 
concentration (l7t l30l|. 

Discussion. The picture of the DNA overstretching 
transition is getting richer and richer as the number of 
experiments increases. However, as we stated in the In- 
troduction, the microscopic nature of the transition is 
still controversial. An interesting approach to the prob- 
lem is to study single-stranded DNA [3l| . Although the 
elongation features in that case could give some insight 
in explaining the rich elastic properties of DNA, the pres- 
ence of two plateaus in single-stranded DNA found in [3l| 
shows that the elasticity of DNA is much more complex 
than previously thought. It seems difficult to put to- 
gether the features of a double-stranded DNA with those 
of single-stranded DNA, where the force-elongation char- 
acteristics depend on the nature of the specific piece of 
DNA object of elongation. 

Possibly related to this, other experimental works sug- 
gest that DNA melting could be responsible for the 
DNA overstretching features [HHH [Hi, Hi ■ Although 
strong evidence in that direction appears reasonable, 
many doubt still remain in this respect [23M25I l32j . 

The model presented here does not consider any mi- 
croscopic aspects of the transition but rather adopts, in 
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the GL spirit, a mcsoscopic point of view, valid, in prin- 
ciple, for any of the above descriptions. We consider 
the transition between two states (the B-DNA and the 
overstretched state) using a stochastic equation of motion 
with an asymmetric double- well potential in the presence 
of one free parameter connected to biophysical proper- 
ties of the DNA chain and/or of the surrounding sol- 
vent. The model is able to reproduce the equilibrium 
(over)stretching curve with good precision. Moreover, 
the approach used here allows us to study the dynam- 
ics of the transition. It can be used to carry out real- 
istic calculations in DNA dynamics, as in translocation 
under strong forces caused by molecular and nanotech- 
nological motors, and represents a valid tool to perform 
out-of-equilibrium modeling. In this sense the model can 



be expanded in order to include other variables able to 
describe more specific microscopic pictures, such as, for 
example, the twist-stretch coupling previously discussed, 
or melting properties. 
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